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Summary 

We examine the multi-point velocity field for non-Gaussian models as a probe of non-Gaussian 
behavior. The two-point velocity correlation is not a useful indicator of a non-Gaussian density 
field, since it depends only on the power spectrum, even for non-Gaussian models. However, 
we show that the distribution of velocity differences vi — V2, where vi and V2 are measured at 
the points ri and r 2 , respectively, is a good probe of non-Gaussian behavior, in that p(vi — v 2 ) 
tends to be non-Gaussian whenever the density field is non-Gaussian. As an example, we 
examine the behavior of p(vi — v 2 ) for non-Gaussian seed models, in which the density field is 
the convolution of a distribution of points with a set of density profiles. We apply these results 
to the global texture model. 

Key Words: Cosmology - dark matter - galaxies: clustering - galaxies - large-scale Structure 
of the Universe. 
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1 Introduction 



A fundamental problem in modern cosmology is to understand how the large scale structures 
in the observed Universe arose from the primordial energy density fluctuations. The simplest 
hypothesis is that an inflationary de Sitter phase produced Gaussian, scale-invariant adiabatic 
fluctuations leading to the generation of large scale clustered structure grew due to their grav- 
itational instability. This scenario, in a Universe dominated by cold dark matter (CDM), has 
been considered the most succesful theory for the formation of galaxy or cluster structures in 
the Universe (Frenk et al. 1988; Davis et al. 1992). 

However, there are also notable setbacks for this model from observational evidence. The 
most serious of these observational challenges have come from the variance in the cell counts 
of IRAS galaxies (Efstathiou et al. 1990), from the measurement of the skewness in the 
distribution of QDOT-J RAS galaxies out to 140 /i _1 Mpc (Saunders et al. 1991), from the 
large-scale clustering of radio galaxies (Peacock & Nicholson 1991) and rich clusters of galaxies 
(Efstathiou et al. 1992), and from the statistical analysis of the APM galaxy distribution on 
scales beyond 10/i _1 Mpc (Maddox et al. 1990). 

Occurrence of high-redshift quasars (Warren et al. 1987), radio sources (Lilly 1988; Cham- 
bers, Miley & van Breugel 1990) and protoclusters of galaxies (Uson, Bagri & Cornwell 1991), 
as well as the Great Bootes Void, Great Wall, Great and Giant Attractor regions (Kirshner et 
al. 1981; Geller & Huchra 1990; Lynden-Bell et al. 1988; Scaramella et al. 1989) may pose 
further problems to the "standard" CDM Gaussian model (see the discussion in Davis et al. 
1992). 

As a consequence, alternative theories of the origin of large scale structure in the Universe 
continue to be worth investigation. Among the chief rivals to the inflationary model are a 
class of models in which the initial inhomogeneities are produced by a primordial population of 
seeds, e.g. cosmic string loops (Zel'dovich 1980; Vilenkin 1981), global textures (Gooding et al. 
1992 and references therein), or generic "seeds" (Villumsen et al. 1991; Gratsias et al. 1993), 
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and the final linear matter fluctuations are not described by a Gaussian distribution. Non- 
Gaussian distributed perturbations are the most general statistical framework for computing 
cosmological observables, such as spatial galaxy correlation functions (Matarrese, Lucchin & 
Bonometto 1986; Scherrer & Bertschinger 1991), expected size and frequency of high density 
regions (Catelan, Lucchin & Matarrese 1988), hotspots and coldspots in the cosmic microwave 
background (CMB) distribution on the sky (Coles & Barrow 1987; Kung 1993) and higher 
order temperature correlation functions (Luo & Schramm 1993; Gangui et al. 1994). The 
importance of the signature of the initial distribution of non-Gaussian density fluctuations for 
the final clustering pattern of the cosmic structures has been recently emphasized, in a series 
of N-body simulations, by Messina et al. (1990), Moscardini et al. (1991), Matarrese et al. 
(1991) and Weinberg & Cole (1992). Coles et al. (1993) show that the mass distribution in 
models with initially non-Gaussian fluctuations exhibits systematic departures from Gaussian 
behavior on intermediate to large scales. 

A powerful method to probe the matter distribution for non-Gaussian behavior is to examine 
the large-scale bulk velocity field, which is, in the linear regime (i.e. large scales), directly 
related to the density field (Peebles 1980). The problem has been addressed by Scherrer (1992). 
Scherrer, noting that the RMS velocities are the same function of the power spectrum for 
both Gaussian and non-Gaussian density fields, examined the distribution of velocities for 
non-Gaussian density fields. For density fields in which the gravitational potential field is 
non-Gaussian by construction ("local" non-Gaussian fields) the velocity field can be strongly 
non-Gaussian. However, Scherrer found that the one-point distribution of linear velocities for 
"seed" models, in which the density field is the convolution of a discrete set of points with a set 
of density profiles, is highly Gaussian even when the density field is strongly non-Gaussian; this 
result has also been seen in numerical simulations of the texture model (Gooding, et al. 1992). 
The reason for this is the central limit theorem: the linear velocity field is an integral over the 
density field; if the density field is sufficiently uncorrelated, the sum of uncorrected densities 
produces a Gaussian velocity field. Of course, the linear velocity field in this case cannot be 
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an exact multivariate Gaussian, as the divergence of such a field (which is proportional to the 
density field) would then have to be exactly a Gaussian. 

In this paper, we exploit this fact to derive a velocity statistic which is sensitive to non- 
Gaussian density fields. We extend the analysis of Scherrer (1992) to derive the two-point 
velocity distribution for seed models. We show that the distribution of velocity differences is 
an effective probe of non-Gaussian behavior in the density field. 

In the next section, we review the general properties of linear cosmological velocity fields, 
valid for any mass density distribution, Gaussian or non-Gaussian. In §3, after reviewing the 
statistical properties of seed models, we work out exact analytical expressions for the 1-point 
and 2-point velocity distributions using the partition function for randomly distributed seed- 
objects. These results are used to derive an expression for the distribution of the difference 
between velocities measured at two different points. We apply our results first to a toy seed 
model to illustrate some of our general arguments about the usefulness of this method, and 
then to the global texture model with cold dark matter (Gooding et al. 1992). Our results are 
discussed in Section 4. 



2 Cosmological Velocity Fields 

During the linear regime, when the amplitude of matter fluctuations is very small (5 <C 1), 
the peculiar velocity v is related to the matter fluctuation 5 by the equation (see, e.g., Peebles 
1980) 

V- v(r,t) = -HQ 0S a5{r,t) , (1) 

where H is the Hubble parameter, a is the scale factor, and Q is the mean density in units of 
the critical density. Any transverse mode, defined by the condition V • v = and corresponding 
to the rotational part of the velocity field, decays with the expansion (namely, the vorticity 
decays with rate oc a -2 in an Einstein-de Sitter universe). Therefore, the integral expression 
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of the solution of eq. (1) may be written as (Peebles 1980) 




2 n-° A 
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V0(r,t), 



(2) 



which is manifestly irrotational; is the peculiar gravitational potential. We see that in the 
linear regime, both density and velocity fields may be derived from a unique potential <fi . In 
eqs. (l)-(2), we consider only the growing modes for density and velocity. In what follows, we 
fix t equal to the present time, taking a = 1 and H = H . 

Results for the linear velocity field are of interest because the observations, when smoothed 
on sufficiently large scales, recover this linear velocity field. This smoothing is usually done by 
filtering on the scale R the velocity field v(r) by means of a window function W R (x), x = |x|, 



The windowing convolution averages the velocity field v over the point x in a volume ~ R 3 , 
where one has no interest in the substructure. Typical choices are the Gaussian and the top- 
hat window functions (e.g. Bardeen et al. 1986). Working with the Fourier transforms of the 
quantities in eq. (3), it is easy to show that the smoothed linear velocity field corresponding to a 
particular density field is identical to the unsmoothed linear velocity field corresponding to that 
density field smoothed with the same window function (Scherrer 1992). Thus, the smoothed 
velocity field can be derived by applying eq. (2) to the smoothed density field. 

An equivalent manner to describe the velocity field is in terms of the Fourier components 
of the density field, S(k). An important cosmological observable, the rms peculiar velocity a v , 
can be expressed very simply in terms of the power spectrum P{k): 



velocity a v and the power spectrum P{k) is valid for any type of density field, Gaussian or 
non-Gaussian (Scherrer 1992). Thus a v cannot be used to distinguish Gaussian from non- 
Gaussian models. For this reason Scherrer (1992) was motivated to examine the distribution 




(3) 




(4) 



where P{k) is defined by (5(k!)5(k 2 )) = 5o(ki + k 2 )P(/ci). The relation between the rms 
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of v, which is different for Gaussian and non-Gaussian models. For example, in local non- 
Gaussian models, in which the potential field is local nonlinear function of a Gaussian field 
(Kofman et al. 1989; Moscardini et al. 1991) the one-point distribution of a single component 
of v is strongly non-Gaussian (Scherrer 1992). However, for seed models, in which the density 
field is the convolution of a density profile with a random distribution of points, the distribution 
of a single component of v is nearly Gaussian. 

We may look for better tracers of a non-Gaussian density field among the multi-point 
velocity distribution. An obvious choice is the two-point velocity correlation function 



We see that just as in the case of the rms velocity, the two-point velocity correlation depends 
only on the power spectrum. Thus, it is incapable of distinguishing a Gaussian density field 
from a non-Gaussian field having the same power spectrum. 

We argue here that a better test for non-Gaussian behavior is the distribution of velocity 
differences p(vi — V2), where vi and V2 are the velocities measured at the points ri and r 2 , 
respectively. In the usual way, we define v\\ and v± to be the components of the velocity parallel 
and perpendicular to ri — r 2 , respectively. Then it is plausible that the distribution of vi\\ —V2\\ 
will become non-Gaussian for a non-Gaussian density field in the limit where |ri — r 2 | — > 0. [We 
measure vi\\ — v 2 \\ relative to a coordinate system in which the direction of increasing distance 
points from r 2 to 1*1; this insures that vi\\ — i> 2 || has the same sign as dv x /dx in the limit 
|ri — r 2 | — > 0.] From eq. (1), we know that the distributions of V • v and 5 must be the same 
(up to a multiplicative constant). For an isotropic density field, dv x /dx, dv y /dy, and dv z /dz 
must all have the same distribution. There is no simple relationship between the distribution 
of dv x /dx and the distribution of S, since the various partial derivatives of v are correlated. 
However, it would require a very contrived density field for dv x /dx, dv y /dy and dv z /dz to all 
have a nearly Gaussian one-point distribution, while S was highly non-Gaussian. Furthermore, 
vi\\ — V2\\ is proportional to dv x /dx in the limit where the separation between ri and r 2 goes 




(5) 
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to zero. Thus, we expect p(vi\\ — V2\\) to be non-Gaussian for sufficiently small separations 
whenever p(5) is non-Gaussian. This is our motivation for examining p(vi — v 2 ) as a probe of 
non-Gaussian behavior. In essence, the taking of velocity differences is a poor man's derivative. 
Integrating over the density field to derive the velocity field can drive p(v) to a Gaussian because 
of the central limit theorem; taking velocity difference restores the non-Gaussian behavior lost 
through integration. Of course, a major issue is the maximum separation over which non- 
Gaussian behavior can be detected in the distribution of Vi — v 2 ; if p(v) is Gaussian, and if 
velocity correlations go to zero at sufficiently large distances, then the p(vi — v 2 ) must also go 
to a Gaussian over those distances. 



3 Velocity Differences in Seed Models 

To test our velocity-difference statistic, we analyze the large-scale velocity field in a special class 
of intrinsically non-Gaussian models, in which the initial linear density field can be described 
as the convolution of a set of density profiles with a random distribution of points (Scherrer & 
Bertschinger 1991; Scherrer 1992). This set of non-Gaussian models provides a useful test of 
our method for several reasons. Such distributions are analytically tractable and can describe, 
for example, the density field produced by textures (Gooding et al. 1992) and the seeded hot 
dark matter model (Villumsen et al. 1991; Gratsias et al. 1993). Further, such models tend 
to produce a Gaussian one-point velocity distribution even when the density field is strongly 
non-Gaussian (Scherrer 1992), so they are useful to determine if the velocity differences provide 
a better probe of the non-Gaussian behavior. 

In these seed models, the linear density field 5(r) may be written as a convolution of a given 
density profile f(m,r) with a set of points having mean density n(m), where m is simply a 
parameter or set of parameters such as mass or formation time: 

S ( T ) = £ /("ifc, r - x h ) (6) 
h 
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where the sum is taken over all seeds with positions x^. The statistical properties of the density 
field are entirely determined by f(m,r), n(m), and the spatial correlations of the seeds (see 
Scherrer & Bertschinger 1991 for more details). Here we follow Scherrer (1992) and consider 
only randomly distributed seeds for which / is spherically symmetric: /(to, r) = /(to, r). This 
simpler model is a reasonably good description of the texture model (Gooding et al. 1992) and 
the seeded hot dark matter model (Villumsen et al. 1991; Gratsias et al. 1993). 

In the linear regime, both the density perturbations and the velocities induced by the seeds 
can be added linearly, so that the perturbation 5(r) in eq. (1) produces at r the peculiar velocity 
v(r) given by 

vf^^vKr-xJ, (7) 

h 

where v(to, r — x) is the contribution to v, measured at the point r, due to a single m-seed 
located at the point x. From eq. (1), we have 

V • v(m, r - x) = -H fi a6 /(to, |r - x|) ; (8) 

Thus, /(to, |r — x|) acts exactly as a density perturbation source for the single seed-velocity 
contribution v(m, r — x), and the total velocity v(r) is the superposition, for the entire popula- 
tion of seeds, of the single seed velocities. To determine v(m, r — x), we integrate the previous 
equation on the sphere V* centered in x. Gauss's law gives 

— r — x 

v(to, r - x) = -/3 /(to, |r - x|) _ ^ , (9) 

where f3 = H f2 - 6 : we see that the m-seed at x produces a velocity component v(m, r — x) at 
r which points toward the seed location, i.e., each seed tends to create a spherical velocity field 
around it. In eq. (9), we have defined the quantity 

/( m ,| r _x|) = — L— /'"Vr' 2 /Kr'), (10) 
|r — x| z Jo 

which represents the integrated density perturbation due to the seed at x, enclosed in a sphere 
of radius |r — x| centered on the seed. From Eqs.(7) and (9) we obtain the total velocity v in 



term of the accretion pattern, 

v(r) = f(m h , |r - x fc |) ±^*L (11) 

Equations (7)-(ll), are somewhat more general than those given by Scherrer (1992), in that 
the location of the seed is arbitrarily chosen. 

We now analyze the statistical properties of the velocity field v for these seed models. If 
the seeds are randomly distributed, the probability distribution of the total velocity v may be 
obtained by applying the Poisson model (Peebles 1980; Fry 1985). To do this, we introduce the 
stochastic Poisson variables in such a way that 



1 if a seed G cbt/i 

(12) 

otherwise , 



with mean 

(e h ) = n(m h ) dm h dx h . (13) 

From eq.(7), we know that v(r) is the superposition of all the single seed contributions, and to 
remember the stochastic information, i.e., presence or absence of a seed in cbc/i, we define the 
variable v(r) = Eft v(m A , r - x /l )e ft . 

The probability distribution function p(y) may be calculated according to the definition 
(see, e.g., Ma 1985) 

I r+oo 

p(v) = (Mv(r) - v(r)] ) = -— / dH e-"-« LK(t) , (14) 

[Z7T) J-oo h 

where eh (t) is the characteristic (or moment generating) function for the single seed discrete 
process (Peebles 1980; Fry 1985) 

(j) eh (t) = (e-^v^.r-x,)^ = exp j n ^^ dmh £ Xh ^it-v(m h ,r- Xh ) _ ^ 

Therefore, the characteristic function for the total velocity is the product of the individual <p eh 

0(t) = n<Mt) — * ex P / n(m) dmd 3 x ( e lt ^ m ' r - x ) - l) (16) 
h J 
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in the continuum limit drrih — > and rfx^ — > 0. Finally, the velocity pdf p(v) may be written 



as 

p(v) = — l — [ + °° dHe~ itv exp f n{m)dmd 3 x(e it ^ m ' r -^ - l) , (17) 

(271") d J-oo J ^ ' 

which is a generalization of the result given by Scherrer (1992). 

Scherrer (1992) showed that the distribution of the one point velocity, p(v), is a rather 
poor indicator of non-Gaussian behavior for seed models with randomly distributed seeds. In 
particular, for the seeded hot dark matter model, a non-Gaussian velocity field requires an 
extremely low seed density, while the texture model with cold dark matter has a velocity field 
which is nearly indistinguishable from a Gaussian. 

We now consider the two-point distribution of velocities, p(vi,v 2 ). Taking advantage of 
our previous results, it is straightforward to show that the the 2-point distribution is given by 

p(vi, v 2 ) = ( <$r,[v(ri) - v(ri)] 5 D [v(r 2 ) - v(r 2 )] ) 

-— 6 d% dh 2 e-ti-v!-it 2 .v 2 0(ti? ts) ? (lg) 



(27T) 

where 

0(ti,t 2 ) = exp 



J n(m) dmd 3 x (e itl,v(m ' ri - x)+it2 - v(m>Pa - x) - l) 



(19) 



is the characteristic function for the 2-point distribution of velocities; it generates the 2- 
point velocity correlation tensor component defined by (Groth, Juszkiewicz & Ostriker 1989) 
£a/3(i"i,r 2 ) = ( v a (ri) vp(r 2 ) ) , a, (5 — 1, 2, 3 , through partial differentiation: 

W(ri,r 2 ) = -|^(0)|^(0). (20) 

The statistical properties of the large-scale peculiar velocity field, as described by the correlation 
tensor, are extensively examined in Gorsky et al. (1989) and Landy & Szalay (1992). 

From eq. (18), one can in particular obtain the probability distribution of the velocity 
difference 

P(vi - v 2 ) = j^y 3 j ^ dh e-^-^-v.) m (21) 
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where now 

0(t) = exp J n{m) rfmrf 3 x (e lt ^ (rn ' ri - x) -^ (m ' r2 ' x)1 - l) . (22) 

A convenient measure of the deviation from a Gaussian is given by the cumulants /%>, defined 
by 

W(i) = E^ (23) 
P =i P- 

For the distribution of velocity differences, it is clear that 

k p = J n(m) dm d 3 x. [v a (m, ri — x) — v a (m, r 2 — x)] p (24) 

where eq. (24) gives the cumulants corresponding to the distribution of the a component of 
the velocity differences. The normalized cumulants 

X p = k p /k p 2 /2 (25) 

give a measure of the deviation of a distribution from a Gaussian; for \ p <^ 1 the distribution 
is nearly Gaussian, while \ p ^> 1 indicates a highly non-Gaussian distribution. 

We argued in §2 that the distribution of vi — v 2 should show the same level of non-Gaussian 
behavior as the distribution of 5; here we test these claims numerically. First consider a toy 
seed model which consists of randomly-distributed seeds with number density n , all having 
the same mass m s , so that n{m) = n 5ri{m — m s ). We also assume a spherical tophat density 
profile: f(m,r) = mg for r < r , and /(r) = for r > r . [This is not a physically realistic 
model, but the distributions for 5 and v can be derived analytically, allowing us to confirm 
some of our general arguments about the usefulness of the velocity difference statistic. Some 
of the properties of this model were examined by Scherrer & Bertschinger (1991). Alternately, 
this model gives the velocity field smoothed with a spherical tophat window function for a 
set of point-like seed masses]. To illustrate the points discussed earlier, we calculate \ p for 
the distributions of 5, v x (a single component of the velocity) and dv x /dx. If we define n = 
n 47rrQ/3, then the distribution of S is just a Poisson distribution with 

A p = n l ' p/2 (26) 
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The value of A p for v x can be derived using the techniques outlined in Scherrer (1992); we find 
(for even p) 

^ = «'^(?TT)(?T3 + *b)(§r (27) 

Finally, consider the distribution of dv x /dx, which is proportional to — v 2 \\ when the sepa- 
ration goes to 0. In this case, eq. (24) with v a (m, r 1 — x) — v a (m, r 2 — x) replaced by dv x /dx 
gives 

A p = n^iB/Qy^-iy + f^x 2 - lfdx\. (28) 

The lowest order non-zero deviation from a Gaussian for v x is given by p = 4. For this case 
we find that A 4 = n" 1 (distribution of 5), A 4 = 0.14n _1 (distribution of v x ), and A 4 = 0.45n _1 
(distribution of dv x /dx). Furthermore, in the limit where p — > oo, the asymptotic expressions 
for A p are 

^-^({f- (29) 

for v x , and 

x ^- vr w^T' <3o) 

for dv x /dx. Thus, for an appropriate choice of n, the distribution of v x can be nearly Gaussian 
while the distribution of S is highly non-Gaussian, in agreement with the results of Scherrer 
(1992). However, the distribution of dv x / dx displays roughly the same deviation from Gaussian- 
ity as does the distribution of 5. Therefore, for sufficiently small separations, the distribution 
of — i> 2 || will also trace the non-Gaussian behavior of 5. 

Now consider a physically realistic non-Gaussian model: the Q = 1 CDM global texture 
model (Gooding et al. 1992 and references therein ). This model can be approximated as a seed 
model, with a random distribution of seeds (the locations of the texture knots). We take the 
parameter m to be the conformal time r at which the knots unwind; then the number density 
of knots unwinding at r is (Gooding et al. 1992) 

dn = uT~ 4 dr, (31) 
13 



with v = 0.04. We use an analytic approximation for /(r, r) given by Gooding (1992) 

SM = C °-^- r ^ S „ (32) 

where A = 1.88 and the constant C does not affect our calculations. The growth factor f g is 

f g = a^ia^S^a) - 5 1 (a i )5 2 (a)), (33) 

where a is the scale factor at which the density perturbations are measured, is the scale factor 
corresponding to the conformal time r, and Si and <5 2 are the growing and decaying modes in 
linear theory (see Scherrer 1992 for more details). 

Using this approximation for the CDM texture model, we have numerically evaluated the 
normalized skewness A3 of the distribution of the component of the velocity difference along 
the separation ri — r 2 , v\\(ri) — v\\(r 2 ) , as a function of separation distance |ri — r 2 |. In Figure 
1, we plot this skewness versus the separation distance of the points at which the velocities 
are measured (solid curve). We see that for relatively small separations, the difference of 
velocities is a good tracer of non-Gaussian behavior. Our statistic is useful up to a separation 
of about I0h~ 2 Mpc. The density field for the texture model has positive skewness, so it is 
not surprising that that i>||(ri) — W||(r 2 ) has negative skewness (see eq. 1 and our argument in 
§2). It is trivial to extend these results to the case of the smoothed velocity field because, as 
noted in §2, the smoothed linear velocity field is identical to the unsmoothed linear velocity 
field corresponding to the smoothed density field. For a velocity field smoothed with a Gaussian 
window function W(r) = exp(— r 2 /rg), the skewness of the distribution of velocity differences 
along the separation vector is given in Figure 1 for r = 2hr 2 Mpc (dashed curve). We find, 
surprisingly, that the smoothed velocity field actually deviates more strongly from a Gaussian 
than the unsmoothed field. The density fluctuation produced by a single texture integrated out 
to a radius R goes to zero in the limit R ^> r (Gooding, et al. 1992) so smoothing a texture 
with a window function of radius r 3> r effectively erases the texture, and we are left with a 
density field consisting of larger textures with larger separations, giving a more non-Gaussian 
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field. Thus, the enhancement of the non-Gaussian nature of p(v\\(ri) — v\\(r 2 )) with smoothing 
is probably peculiar to the texture model. 

4 Discussion 

Our results indicate that the distribution of f||(ri) — i>||(r 2 ), the component of the velocity 
difference at two points along the vector separating the observation points, is a useful probe of 
non-Gaussian behavior. This distribution tends to be non-Gaussian whenever the underlying 
density field is non-Gaussian, in contrast to the one-point distribution of velocities, which 
tends to be more Gaussian than the underlying density field (Scherrer 1992). Although we have 
examined only a single class of non-Gaussian models, our results should be generally applicable 
to any non-Gaussian model, because of our argument in §2. We note that a related argument 
for detecting signatures of non-Gaussian density fields in the cosmic microwave background 
has been given by Moessner, Perivolaropoulos, and Brandenberger (1994). They argue that the 
temperature gradient of the CMB should provide a better test to distinguish the cosmic string 
model then the distribution of the temperature field itself. 

P.C. is grateful to S. Matarrese for helpful discussions. P.C. was supported in part by the 
Italian Ministero della Ricerca Scientifica e Tecnologica and by the Fondazione Angelo Delia 
Riccia. R.J.S. was supported in part by the DOE (DE-AC02-76ER01545). 
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Figure Captions 



Figure 1. The normalized skewness A 3 of the distribution of f||(ri) — f||(r 2 ) [the component of 
the velocity difference parallel to the separation ri — r 2 ], as a function of the separation |ri — r 2 | 
for the unsmoothed density field (solid curve) and the density field smoothed with a Gaussian 
window function of radius 2h~ 2 Mpc (dashed curve). 
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